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Abstract. The MIMAC project aims to provide a nominal fluorine de- 
tector for directional detection of galactic dark matter recoil events. 
" Its expected behavior reaches an important part of the predicted spin 

dependent elastic scattering interactions of the supersymmetric neu- 
tralino with protons. Hence, the parameter space in the MSSM and 
the NMSSM models with neutralino dark matter could be probed by 
J^. such experimental efforts. In particular, a good sensitivity to spin 

^ dependent interactions tackles parameter space regions to which the 

predictions on spin independent interactions and indirect signatures 
I— I are far below current and projected experiments. 

> 

1 Introduction 

\^ Recent and rapid experimental developments have been constraining particle dark 
matter (DM) candidates. In particular, the direct detection community has showed 

^\ its capacity to scan DM-nucleon interactions for candidates of masses at the elec- 
troweak scale with spin independent interactions down to ~ lO"^** cm^ and spin 
dependent cross sections down to ^ lO^''^ cm'^ (see (Censier, 2011 1 for a review 
[ on the subject). 

^ The directional technique relies on the fact that the solar system is in motion with 
respect to the galactic reference frame with a velocity pointing toward the Cygnus 
constellation. The interaction between the DM of the galaxy and a detector on 
Earth would happen in a preferred direction: recoil events could record this asym- 



X 



metry. In (Billard et al., 2010a Billard et al., 2010b) the projected sensitivity 



to spin dependent proton-DM interactions of a forthcoming fluorine detector is 
estimated. The characteristics of the simulated detector are set to be those of the 
MIMAC project, but are also representative of the whole generation of detectors 
currently in development. Data was simulated and analyzed for a detector made 
of 10 kg of CF4, operated at 50 mbar, assuming that recoil tracks can be solved. 
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including the head-tail determination of the event, with a threshold at 5 kcV, and 
an exposure of 30 kg yr. 

Supersymmetric models provide neutralinos as a DM candidate. The Minimal Su- 



persymmetric Standard Model (MSSM) and Next-to-MSSM (NMSSM) (Ellwanger 



et al. 



est supersymmetric particle (LSP) is stable, 
non-universal gaugino masses one could find neutralinos as light as 



2010 ) are well established realizations of Supersymmetry in which the light- 
In the MSSM and NMSSM with 

13 GeV 

(see (Albornoz Vasquez et al., 2011b I and references therein) and ~ 1 GeV 
(see (Albornoz Vasquez et al., 2010 Draper et al., 2011) and references therein) 
respectively, and still overcome all particle physics constraints. That range being 
not far from the lower limit of the direct and directional detection sensitivity, it 
is crucial to make a thorough exploration of the different possibilities, and in par- 
ticular, to establish the range of elastic scattering and annihilation interactions to 
be expected. 



2 Directional detection sensitivity to supersymmetric configurations 

2.1 Exploration of supersymmetric configurations and constraints 

Supersymmetry offers a large number of possibilities in a multidimensional pa- 
rameter space with new interactions and particles. The phenomenology related 
to a neutralino DM candidate is mostly determined by its mass and composition, 
as well as by the Higgs sector. Hence among the relevant free parameters are 
gaugino masses Mi, M2 and M3, the ratio of the vacuum expectation values of 
the two Higgs doublets tan/3, and the fi mass term. The MSSM and the NMSSM 
differ in the Higgs sector, where the NMSSM has an extra singlet scalar field. This 
gives a natural explanation to the energy scale of the MSSM /i mass term, which 
becomes effective in the NMSSM. To complete the Higgs sector parametrization, 
in the MSSM we have the pseudoscalar mass Ma which is replaced in the NMSSM 
by the scalar couplings A and k, as well as the corresponding trilinear couplings 
Aa and A„. The Higgs particle spectrum is expanded from the MSSM -with h, H, 
H=^ and A- to the NMSSM -with Hi, H2, H3, H^^, Ai and A2. In particular, the 
stringent constraints on the lightest MSSM scalar Higgs h do not necessarily apply 
to the lightest NMSSM scalar Higgs Hi, since it could be strongly dominated by 
a singlet component. This, in turn, is at the origin of a broader set of possibilities 
for physics below the 100 GeV scale in the NMSSM. The set of free parameters is 
completed by the soft sfermion masses and the trilinear coupling of the top sector 
At -the other trilinear couplings being set to zero. 

Scanning the multidimensional supersymmetric parameter spaces is not an easy 
task. To this end a Markov Chain Monte-Carlo code was developed to scan 



the multidimensional supersymmetric parameter spaces (Albornoz Vasquez et al.. 



2010) . This code, built on micrOMEGAs 2.4 ( [Belanger et al., 2005[ ), evaluates each 



supersymmetric parametrization using a likelihood function to fit particle physics 
limits, including limits on masses -such as the chargino or the Higgs bosons-, elec- 
troweak observables -such as (g — 2)^ or the Z invisible width-, and B-physics 
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-such as the b — ?> S7, Bg — > /i^/i" and B ^ rvr branching ratios. NeutraHnos 
are also required to represent at least 10% of the cosmological DM as measured 
by WMAP (Komatsu et al., 2011 1, hence a likelihood function was established for 
the relic density of the neutralinos after thermal freeze-out. An iteration process 
based on a Metropolis-Hastings algorithm generates a random walk in the multi- 
dimensional parameter space. 

Once the parameter space is explored and configurations allowed by particle physics 
experiments are found, we may compare the direct and indirect detection yield. 
Therefore, we compute the spin independent cross sections times neutralino-to- 
DM fraction to which we apply the exclusion limits by XENONIOO (Aprile 



et al., 2011 ). We also apply Fermi-LAT constraints on 7-ray fluxes from the Draco 
dwarf spheroidal galaxy (dSph) (Abdo et al., 2010). The 7-ray flux stemming 
from neutralino annihilations in the DM dominated Draco galaxy was estimated 
by computing the 7-ray production cross section and spectrum from neutralino 
annihilations, and multiplied by the line-of-sight integral, provided by Fermi-LAT 
for a Navarro- Frenk- White halo profile (Navarro et al., 1996). For more details, 
see (Albornoz Vasquez et al., 2011a). 



2.2 The impact of directional detection on neutralinos 

Directional detectors based on fluorine are sensitive mostly to the spin depen- 
dent interactions of DM with protons. The directional detection prospects for a 
MIMAC-like detector show that if no background events are recorded, the detec- 
tor would be able to exclude cross section down to ~ 4 x lO"**^ cm^ for a mass 
of 10 GeV. Furthermore, determining the directionality increases the number of 
observables. Using a likelihood analysis, it is possible to determine the WIMP 
mass and interaction cross section. The needed statistics to solve these charac- 
teristics can be translated into a sensitivity curve, which is, of course, above the 
exclusion limit. In light of these sensitivity curves, three regions are defined in the 
^fjp^ vs. rriyO^ plane being the neutralino to DM fraction at Earth's position): 
above the discovery limit, between discovery and exclusion limits, and below the 
exclusion limits. In the first region, neutralinos are expected to be detected, their 
mass and cross section could be measured. Of course, if no event is measured, 
then the corresponding configuration would be ruled out. In the second region, 
neutralinos would produce some signal but not enough to be solved. If they do 
not, then scenarios lying in this region would be excluded. Finally the third region 
is populated by neutralinos that would not yield any effect in the forthcoming di- 
rectional detection experiments. 

We present the results for the MSSM and NMSSM in Fig. [l] It is important 
to notice that many configurations lie above the projected exclusion limit of MI- 
MAC in both the MSSM and the NMSSM. This already encourages the efforts for 
building fluorine directional detectors. We show in green points allowed by both 
XENONIOO and Fermi-LAT, while points failing one or the other are tagged in 
yellow. Those excluded by both are tagged in red. 

The predictions for neutralinos in these two models only differ significantly for neu- 
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Fig. 1. Spin dependent proton-neutralino interactions as a function of neutralino mass 
in the MSSM (top panel) and the NMSSM (bottom panel). The projected sensitivity 
of a typical directional detector are also shown. In green: points safe with respect to 
XENONIOO limits on spin independent interactions and Fermi-LAT limits on 7-rays from 
the Draco dwarf spheroidal galaxy. In yellow: points in conflict with either one or the 
other. In red: points failing to overcome both exclusion limits. 



tralino masses below ^ 30 GeV, This is a direct consequence of the very different 
configurations in the Higgs sector 



Indeed, as it was found in (Albornoz Vasquez 
et al., 2010] ) and further explored in ( [Albornoz Vasquez et al., 2011a| ), NMSSM 



neutralinos below 30 GeV achieve the relic density by resonantly annihilating via 
light scalar or pseudoscalar Higgs bosons, absent in the MSSM. In the latter, the 
only possibilities are to push as much as possible the masses of the lightest scalar 
Higgs in order to allow lighter pseudoscalar masses, and have efficient enough anni- 
hilations. This scenario is actually heavily constrained and should not be regarded 
as plausible, as shown in (Albornoz Vasquez et al., 2011b). The other possibility is 
exchanging very light staus, which provides neutralinos down to 12.6 GeV. Since 
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the MSSM is contained in the NMSSM, this configuration is also realizable in the 
NMSSM. Hence, a detection below 30 GeV implies different predictions for both 
models: in the MSSM we would expect to observe light sleptons (< 100 GeV), 
while the NMSSM could predict a very light Higgs boson. 

In both models, there is a prominent concentration of points with m^o ~45 GeV, 
corresponding to neutralinos in resonant annihilations through a Z boson, an ef- 
ficient mechanism to attain the correct relic density. The points that can be 
detected around the Z resonance are those not falling exactly in it, but rather 
those of masses < 40 GeV and > 50 GeV. This is easily understood in terms of 
the elastic scattering cross section: the spin dependent interactions occur mainly 
via the exchange of a Z boson. While the Z resonance represents a good way to 
obtain a plausible relic density, a too fine-tuned relation between neutralino and Z 
masses leads to too small abundances unless the coupling to the Z is small. Thus, 
for points with m^o sitting too close to Mz/2, the ZxiXi coupling is small, hence 
spin dependent interactions are consequently reduced. 

For the larger masses, for LSP with larger higgsino components, the "ZiXiXi cou- 
pling is usually large, since it is proportional to A'';^3 — N^^ {Nf^ and N^^ being 
the higgsino-d and higgsino-u fractions of the neutralino). However, above a few 
hundred GeV the mass split between the lightest neutralino and squarks narrows. 
Hence both Z and g^, d contribute to the interactions. It turns out that these two 
contributions are destructive. Therefore, for a generally dominating Z exchange 
with rather large couplings, those configurations having a large enough squark 
exchange can lower the spin dependent proton-neutralino cross section by a few 
orders of magnitude. This is why not all the higgsino points have large interac- 
tions. Thus, only a fraction of them falls in the discovery region. The general trend 
to have smaller cross sections towards larger neutralino masses is a consequence of 
the kinematic behavior of the cross section: when m^a ^ nip, then the neutralino- 
proton cross section is proportional to m~o . When the maximum Z^iXi coupling 
is achieved -i.e., TVfg — Nf^ ~ 0.5-, the upper limit for the interaction cross section 
as a function of the neutralino mass is the -2 power law observed in both panels 
in Fig.[l] 

Hence, in a very general way, any detection of a neutralino implies a large higgsino 
fraction, for a neutralino mass below ^ 150 GeV. This implies that the detection 
of a neutralino fixes the fj, mass term to be lighter than 200 GeV. In turn, this 
predicts a chargino of mass below 200 GeV, which should be observed at LHC. 
The other possibility would be small squark masses, which, however, seem to be 
not an option any more after the first round of results at ATLAS and CMS. 
It is also important to notice that many points escape both the discovery region 
and the exclusion region. This should not be taken as a drawback, but encourage 
experimental efforts to develop as good techniques as possible. Consequently, it 
is also important to keep an eye on other detection techniques, in order to tackle 
the most difficult configurations. 
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2.3 Complementarity of direct, indirect and directional detection techniques 

In Figs. [2] we display the points in the 7-ray flux vs. £_a^^ plane, using the ma- 
genta, cyan and blue tagging for those configurations that could be discovered, 
excluded or would not yield any effect in a nominal directional detector. 
The XENONIOO exclusion limits rule out part of the parameter space, which is 
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Fig. 2. Predicted 7-ray fluxes from the Draco dwarf galaxy as a function of ^cr in the 
MSSM (top panel) and the NMSSM (bottom panel). Points excluded by XENONIOO 
are not drawn. The Fermi-LAT limits for the flux are shown. In magenta: points falling 
in the MIMAC discovery region. In cyan: points falling in the MIMAC exclusion region. 
In blue: points beyond the MIMAC sensitivity. 

not shown in these figures. Hence, a fiuorine based directional detector would scan 
configurations that are safe with respect to XENONIOO limits on spin independent 
interactions. Furthermore, in the MSSM, the magenta point with the smallest spin 
dependent interaction reads £,a^^ ~ 4 x 10"*^ cm'^, which is unlikely to be within 
the sensitivity of the projected XENONlton or other spin independent-oriented 
projected detectors. In the NMSSM this is even more drastic: there is a point in 



MIMAC and neutralinos 



7 



magenta with C""^^ — 10~^° cm^! 

Regarding the indirect 7-ray flux from the Draco dSph, the conclusion is similar: 
we find discoverable configurations which lie up to four orders of magnitude below 
the Fermi-LAT Hmits in the MSSM, and even more in the NMSSM. 
It is important for the prospect of directional experiments that we find large con- 
centrations of points which are not excluded by any experiment yet, which are far 
away from detectability by other techniques such as indirect detection and direct 
detection, and which could be discovered or excluded by such projected detectors. 

3 Conclusions 

The projected sensitivity for directional detectors such as MIMAC would allow to 

probe a large portion of parameter space of neutralino DM supersymmetric config- 
urations, especially towards the lightest LSP regions, below 30 GeV. A detection 
could happen below 150 GeV, and would imply a significant higgsino fraction in 
the neutralino composition, which in turn predicts a chargino lighter than 200 
GeV. For discoveries of even lighter neutralinos, the predictions of the MSSM and 
the NMSSM are quite different: the former points towards light sleptons while the 
latter implies light scalar or pseiidoscalar Higgs bosons. 

The interplay between the projected sensitivity of fluorine directional detectors, 
the spin independent interactions and the 7-ray fluxes expected for neutralinos 
in the MSSM and the NMSSM is a crucial feature for the future explorations of 
neutralino DM. If the LHC tells us something about Supersymmetry, then we may 
have indications for which technique is the most adapted to discover or exclude 
the existence of a neutralino in galactic systems. Conversely, signals in direct, 
directional or indirect detection could help the LHC to confirm or rule out the 
MSSM and/or the NMSSM. 
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